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Abstract 

A study of transitions between T states with the emission of charged pions using 398 fb“^ of data 
collected with the Belle detector at the KEKB asymmetric energy e'’'e“ collider is presented. A 
clear peak from the decay T(15) ^ is observed in the invariant mass distribution of 

pairs from the iJ,~ tt~X) event sample. The mass difference distribution — 

M^+^-) for M^+^- from the T(IS') mass region has two peaks from T(25, 35) ^ T(15)7r'''7r“ 
decays, with no background. A third peak at AM = (1119.0 ± 1.4) MeV/c^ can be interpreted 
as evidence of a signal from the decay T(45) ^ T(15)7r'''7r“ with a subsequent T(15) ^ pC 
transition. This is the first example of a non-BB decay of the T(45) resonance. The preliminary 
estimated branching fraction is equal to .B(T(45) ^ T(15)7r+7r“) = (1.1 ± 0.2(stat.) ±0.4(sys.)) x 
10 - 4 . 

PACS numbers: 
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INTRODUCTION 


The bottomonium state T(4S') has a mass above the threshold for 55 pair production 
and decays mainly into these 5-meson pairs (5(T(4S') —> 55) >96% |l|). The decay 
modes T(4S') ^ T(mS')7r7r with m =1, 2, 3 should exist. These decays are analogous to 
the decay modes of the low lying T states j^. Upper limits on the branching fractions 
of T(4S') ^ T(lS')7r+7r“ and T(4S') —>• T(2S')7r’''7r“ decays have been set by the CLEO 
experiment [^. Similar decay modes of the charmonium state V’(3770), which has a mass 
above the DD production threshold have been observed recently 

In this paper we report the first evidence for the decay mode T(4S') —>■ T(lS')7r’''7r“ from 
the Belle experiment. 


EVENT SELECTION 

In this study 398 fb“^ of data collected by the Belle detector on the T(4S') resonance 
and in the nearby continuum is used. Well reconstructed charged particles and photons 
are used to reconstruct the decay T(4S') —T(lS')7r+7r“ with the subsequent leptonic decay 
T(l^) 

Several selection criteria for charged tracks and neutral particles in an event were applied. 

Charged particle candidates were required to satisfy the following requirements: 

• transverse momentum, pt>30 MeV/c; 

• impact parameter transverse to the beam, dr <3 cm; 

• impact parameter along the beam axis, dz <4 cm; 

• tracks that are not identihed as a decay product of a reconstructed secondary V°. 

Charged particles are identified by combining responses from the CDC, TOF and ACC 
subdetectors of Belle into a likelihood Ci j^, Q where i indicates the particle type (e, fi, 
71, K, p). A charged particle is identihed as an electron if the corresponding likelihood ratio, 
Pe >0.9 j^, or if the electron mass hypothesis has the highest probability. The electron 
detection efficiency for Pg >0.9 is approximately 90% for single electrons embedded onto 
hadronic events and uniformly distributed over the polar angle range 35° < 9 < 125° and 
the momentum range 0.5 GeV/c < p <3.0 GeV/c. Gharged particles are identihed as muons 
if the corresponding muon likelihood ratio is 5^ >0.8. The muon detection efficiency above 
a given likelihood threshold is approximately 91.5% for the single-track simulated muons 
uniformly distributed over the polar angle range 20° < 9 < 155° and the momentum range 
0.7 GeV/c < p <3.0 GeV/c. Gharged particles that are not identihed as muons or have a 
likelihood ratio 5e<0.05 were considered as pions. 

The identihcation of y’s and tt^’s is based on information from the electromagnetic 
calorimeter. Galorimeter clusters not associated with reconstructed charged tracks are con¬ 
sidered as 7 candidates. In addition, these electromagnetic clusters must have energies 
greater than 50 MeV and not satisfy the vr® hypothesis when combined with another photon 
in the event. 

Two electromagnetic showers with an invariant mass — m^o| <8 MeV/c^ (where 
m^o is the nominal mass) form a 7r° if the conhdence level for their kinematic ht to 
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the 7 r° hypothesis is greater than 0.1 and if their reconstructed momentum is greater than 
100 MeV/c. 

Hadronic events are selected by the standard Belle hadronic selection. The most impor¬ 
tant of these selection criteria are the following: multiplicity of charged tracks in an event 
^ch ^3; the electromagnetic calorimeter cluster multiplicity nECL >1; the event’s visible 
energy > 0.2^/s, where ^/s is a center of mass energy; the energy sum of good clus¬ 
ter energies in the electromagnetic calorimeter must satisfy 0.18< < 0 . 8 ; the sum 

of the 2 ; components of each charged track’s and photon’s momenta is required to satisfy 
\P^\<0.5^/s. 

Events with a pair of oppositely charged particles and with an invariant mass 
Mch+ch- >9 GeV/c^ were selected from the hadronic sample. The invariant mass distri¬ 
bution for these events in the on-resonance sample is shown in Fig. 1. The M^+^- invariant 
mass distribution when both particles are muons is also shown in Fig. 1. An enhancement 
near the T(IS') mass can be observed. 

The peak at the invariant mass ~10.6 GeV/c^ originates from the process e^e“ —>■ 
/r'''/r“(e+e“( 7 )). About 60% of the events in the spectrum come from identihed 

pairs. The absence of e+e“ events with high invariant mass is mainly due to selec¬ 
tion criterion for the energy sum of good cluster energies in the electromagnetic calorimeter 
A'sum/\/s<0.8. Therefore we restrict our attention to pairs. 

To reduce the background from poorly reconstructed events we impose an additional 
requirement on the visible energy of selected events: 10.5 GeV <Fvis<12.5 GeV. After all 
cuts described above about 112.5k events were selected. 


ANALYSIS 

To reconstruct T(4S') —> T(lS') 7 r’'' 7 r“ decays, which in addition to the decay products of 
the T(IS') contain two charged pions in the hnal state, events were considered. 

This additional selection criterion reduces the primary sample of events considerably (by 
about a factor of 100). The number of selected events is 957. The muon pair invariant mass 
distribution, for the events is shown in Figs. 1 and 2. The clear low 

background signal for T(IS') —> is shown in Fig. 2. 

The £t to this distribution using a sum of the Grystal Ball function for the signal and a 
polynomial function for the background results in the peak position M^+^- (peak) = (9454.7± 
5.5) MeV/c^, which is consistent with the nominal T(IS') mass value Mx(i 5 ) = (9460.30± 
0.26) MeV/c 2 Q. 

To observe resonance states that decay into the T(IS') tt+tt" hnal state the distribution 
of the mass difference AM = — M^+^-) where M^+^- is restricted to |M^+^- — 

Mr[is) \ <6 MeV/c^ was examined (see Fig. 3). Here My(is) is the nominal TGlS") mass. 
To reduce background from T(IS') production in radiative return processes j9| with the 
subsequent conversion of the emitted photon as an electron-positron pair ( 7 * —>■ e’''e“), 
which is then misidentihed as 7 r’'' 7 r“, we impose an additional requirement on the angle 
between the pion momenta in the lab system, cos^jr-n- <0.95, on the events shown in Fig. 3. 
The cos 6 *^ 7 r distributions for the events in the peaks and in the sideband regions are shown 
in Fig. 4. The peaks near cos^t^tt = 1 are from the radiative return T(IS') background 
described above. The additional requirement removes this background. 

Three peaks are seen in the AM distribution (Fig. 3). The hrst (second, third) peaks 
have values of AM ~ 0.56(0.89,1.12) GeV/c^, respectively. 
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The first and second peaks have little or no background. They originate from the de¬ 
cays T(2S') —> T(lS') 7 r+ 7 r“ and T(3S') —>• T(lS') 7 r+ 7 r“ with a subsequent T(IS') —>• ii~ 
transition, respectively. Fits to the first two peaks using a Gaussian (see Fig. 5a,b) re¬ 
sult in the following values of the peak positions AM(lst) = (562.0 ± 0.1) MeV/c^, 

AM(2nd) = (893.5 ± 0 . 2 ) MeV/c^. The mass difference between T(2S'), T(3S') and 
the T(l^) states presented in PDG [l| is (Mx( 2 S) — Tfx(is)) = (563.0 ± 0.4) MeV/c^, 
(Mx( 35 ) — Mx(i 5 )) = (894.9 ± 0.6) MeV/c^, respectively. The values obtained from the 
hts to the distributions in Fig. 5a,b by Gaussian functions are compatible with the cor¬ 
responding PDG values. The y^/NDF for the hrst and second peaks are 1.4 and 1.8, re¬ 
spectively. The small differences between the mean values obtained from the Gaussian hts 
and the PDG values can be explained by systematics due to the imperfect modeling of the 
AM distributions. We conclude that the hrst and second peaks are produced by the decays 
T(2^) ^ T(15)7r+7r- and T(35) T(l^) 7 r+ 7 r-. 


STUDY OF THE THIRD PEAK 


In contrast to the hrst two peaks, the third peak has a considerably larger back¬ 
ground. The position of the peak is derived from a ht to the distribution in Fig. 6 
using a Gaussian for the signal and a polynomial function for the background to be, 
AM = (1119.0 ± 1.4) MeV/c^, which is in good agreement with the mass diherence 
(Mx( 45 ) — Mx(i 5 )) = (1120.0±3.5) MeV/c^ from the PDG Q. The signal above background 
is determined from the ht to be, iVev = (38.0 ±6.9), with a statistical signihcance of 7.3 stan¬ 
dard deviations. This peak is interpreted as a signal from the decay T(4S') —T(lS') 7 r+ 7 r“ 
with a subsequent T(IS') —transition. 

Using the oh-resonance sample (A/i=10.52 MeV, integrated luminosity / Cdt ~ 40 fb“^), 
the mass diherence distribution shown in Fig. 7 shows only two peaks, which are from T(2S'), 
T(3S') decays. The enhancement at AM =1060 MeV/c^ in this distribution is compatible 
with the background level near the third peak in the on-resonance sample (Fig. 3) and 
lower than the luminosity scaled yield of the third peak in the on-resonance data. For the 
renormalized on-resonance sample the total number of events and background in the interval 
1.11 GeV/c^ < AM <1.135 GeV/c^ corresponding to the third peak is = (5.4 ± 0.8), 
= (1.1 ± 0.4) respectively. The number of events in the shifted interval for the off- 
resonance sample is = (2.0 ±1.4). Unfortunately the low statistics of the off-resonance 
sample preclude using it to subtract background in the on-resonance sample. 

However, additional information can be obtained from the study of the 7 r’'' 7 r“ system. 
The invariant mass distribution of the 7 r+ 7 r“-system for events from the observed 

peaks and background is shown in Figs. 8 a-d. These distributions are directly compared 
with model predictions as the efficiency corrections are negligible. The background, which is 
small, is not subtracted from the M^+^- distributions for the three observed peaks (Figs. 8 a- 
c). 


The distribution for the first peak (T(2S') —T(lS') 7 r+ 7 r“ decay) is well described 

by the Yan model 1^, where the hadronic transition between heavy quarkonia is considered 
as a two-step process: the emission of gluons from heavy quarks and subsequent conversion 
of these gluons into light hadrons. This transition can be described in the context of a 
“multipole” expansion scheme where the gluon fields are expanded in a multipole series. 

The T(3S') —T(15')7r'*~7r~ decay distribution (Fig. 8b) on the other hand is described by 
the Moxhay model [ll|. In this model the decay proceeds through coupling to BB^ B*B*, 
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intermediate states • 


This significant coupling to the fiavored sector is motivated by the fact, that the T(3S') lies 
closer to the open BB production threshold than the T(2S'). This decay mechanism causes 
a double peak structure in the M^+^- distribution. The multipole and coupled-channel 
amplitudes interfere in the Moxhay model. 

The above models were already used successfully by the CLEO experiment l^j to describe 
the Mjr+jr- distributions in the T(2S', SS") —T(lS')7r+7r“ decays. The distribution 

for the decay T(4S') —>• T(lS')7r’''7r“ (Fig. 8c) is better described by the Yan model. 

The distributions for all resonance decays (Fig. 8a-c) show an enhancement at high 
masses. In contrast in the distribution for the background (events out of peak ranges, 

i.e. from the mass difference ranges 0.3 GeV/c^ < AM <0.555 GeV/c^, 0.57 GeV/c^ < 
AM <0.885 GeV/c2, 0.905 GeV/c2< AM<1.11 GeV/c^, 1.135 GeV/c^ < AM<1.4 GeV/c^) 
the high mass region is significantly suppressed (see Fig. 8d). This difference in the be¬ 
haviour of the distribution is an additional argument in favour of a resonance inter¬ 

pretation for the third peak. 

The branching fraction for the T(4S') —> T(lS')7r’''7r“ decay was extracted from 
i3(T(4^) ^ T(l^)7r+7r-) = Aev/(Y*f • £ ■ i3(T(l^) ^ /i+/i")). The total number of 
T(4S') in the data sample is = (386 ± 5) x 10®, the nominal branching fraction 
i3(T(lS') —= 2.48%. The EvtGen event generator □ with a matrix element 
taking into account particle spins was used for the Monte Garlo simulation of T(4S') — 
T(lS')7r+7r“ —/i^/i^TT+Tr" events. The Monte Garlo generated events were then passed 
through the detector simulation and reconstruction programs. The efficiency of the event 
reconstruction is then obtained from the Monte Garlo sample to be £ = 0.035. The sys¬ 
tematic uncertainty on the reconstruction efficiency is estimated by comparing the number 
of reconstructed T(2S', 35') events with the number of events that were calculated for the 
radiative return process e+e“ —^ e’''e “7 —> T(2S', 3S')7 in the model [l^. The difference is 
35% and is a conservative estimate of the event reconstruction uncertainty. Another system¬ 
atic uncertainty comes from the poor knowledge of the T(4S') ^ T(lS')7r+7r“ —> 
decay matrix element. This was estimated from a comparison of the previously calculated 
efficiency £ with the efficiency calculated with a phase space matrix element and this sys¬ 
tematic uncertainty results in 8%. 

The preliminary result for the branching fraction is i3(T(4S') —> T(15)7r+7r“) = 
(1.1 ± 0.2(stat.) ± 0.4(syst.)) x 10“"^. The width, T, of the decay T(4S') —> T(15)7r’''7r“ 
is (2.2±1.0) keV, which is comparable with r(T(2S') ^ T(lS')7r^7r“) = (8.1 ± 1.2) keV and 
r(T(3^) ^ T(l^)7r+7r-) = (1.2 ± 0.2) keV. 


CONCLUSIONS 

A study of transitions between T states with the emission of charged pions at Belle 
has been performed. A clear peak from the T(IS') —decay in the invariant mass 
distribution of (p^p“) pairs from the event sample is observed. The mass 

difference distribution — M^+^-) for M^+^- from the T(IS') mass region has 

two peaks from ‘T(2S', 35') —> T(15')7r+7r“ decays with no background. A third peak at 
AM = (1119.Oil.4) MeV/c^ is interpreted as evidence of a signal from the decay T(45') —>■ 
T(15')7r+7r“ with a subsequent T(15') —>■ transition. The preliminary estimated 

branching fraction results in;B(T(45) ^ T(15')7r+7r“) = (1.1 i0.2(stat.) i0.4(sys.)) x 10“'^. 
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FIG. 1; The invariant mass distribution of charged particle pairs for the standard Belle 

hadronic on-resonance sample. 




FIG. 2: The invariant mass distribution of muon pairs in vr+vr X events fitted with 

the sum of a Crystal Ball and polynomial function. 



FIG. 3: The mass difference AM = — M^+^-) distribution where M^+^- lies in 

the T(IS') mass region. Arrows show the positions of the mass differences (Mx( 25 ) — -^t( 15 ))) 
(Mx( 35 ) — Mx(i 5 )) and (Mx( 4 s) — Mx(i 5 )) from PDG P|, respectively. 
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FIG. 4: The cosOttw distributions for events from the three observed peaks and background. 


11 

















































































FIG. 5; The mass difference AM = distribution for the first (a) and second 

(b) peak (|M^+^- — Mx(i 5 )| <6 MeV/c^). The central parts of the peaks are fitted by Gaussian 
functions. 



FIG. 6: The fit to the third peak in the mass difference AM = {M^+^-^+^- — M^+^-) distribution 
(|M^+^- — Mx(i 5 )| <6 MeV/c^) using the sum of a Gaussian and polynomial function. 
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FIG. 7: The mass difference AM = (M^+^-.^+^- — M^+^-) distribution for the off-resonance 
/i'''//“ 7 r+ 7 r“A-sample (|M^+^- — My(i 5)|<6 MeV/c^). 






FIG. 8 : The 7 r'’' 7 r invariant mass distributions for events from the observed peaks and background. 
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